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Beschreibung 

Verfahren zum . Bereitstellen von Fehlerinformationen Otoer In- 
konsistenzen in. einem System von Differentialgleichuhgen .. 

Die Erfindung betrifft ein Verfahren zum Bereitstellen von 
Fehlsrinformatifcnen uber Inkonsistenzen in einem System von 
Dif f erentialgleichungen. ■ ■ ; ' 

Im Stand der .Technik ist 'es bekannt, ein technisches System 
, durch systeme yon Dif f erentialgleichungen zu.beschreiben.. Un- 
ter z^hilfenahme von numerischen Verfahren mit einem Compu- 
tersyitem oder eines ..Anaiogrechners ist es moglich, diese Sy- 
stems zu sirauliere.n, indem fur bestinmte Start- oder Randbe- 
^ dingungen Losungen der Systeme von Dif f erentialgleichungen 
berechnet warden. Beispiele fUr die Be^chreibung technischer 
/ systeme 'sowie typische Simulationsverfahren konnen aus Stan- 
dardwerken der jeweiligen Fachrichtung entnonuaen werden. Aus- 
/ ■ fuhrliche Anleitiingen hierfur sind beispielsweise in G. 
20 Schmidt, ' "Grundlagen der Regelungstechnik" , 2. Auflage-, 

Springer,. Berlin, 1937, in Unbehauen, "Regelungstechnik I", 
6. Auflage, Friedr. Viehweg & Sohn, Braunschweig/Wiesbaden, 
• . .198 9, in' Unbehauen, "Regelungstechnik II'',. "S. Auflage, - 

Friedr. Viehweg & Sohn, Braunschweig/Wiesbaden,- 1989, in Un- 
2 5 behauen, "Regelungstechnik III", . 2 . Auflage, Friedr. .Viehweg ' 
& Sohn, Braunschweig/Wiesbaden, 1986, in E, Pfeiffer, "Ein- 
' ftlhrung in die Dynami.k", B. G. Teubner . Stuttgart , 1989, sowie 

in E. Ziegler (Hrsgb. ), "Teubrier-Taschenbuch der Mathematik" 
■ B. G. Teubner, Stuttgart, 1996 angegeben, Bei der Simulation 
tritt hkufig das Problem auf, dali ein. Verfahren zum numeri- ' 
schen Losen eines Systems von Dif f erentialgleichungsn ab- 
bricht, well da^ zugrundeliegende System von Differential- 
gieichungen singular ist. Aus demgleichen Grund kann es vor- 
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.kommen,. daii ein Analogrechner* das, Systemverhalten nicht vor-- 
hersehen kann. Zudem kommt es haufig vor, dali.beiin Losen der 
Systeme von Dif f erentiaigleichungen Losungen auftreten, die 
nicht plausibelsind - 
5 ' 

In einem solchen' Fall ist es im Stand der Technik bekannt, . 
daB der Prograiniuierer oder ein Benutzer eines Prograrams zum 
Losen der Dif f erentialgleichungen mit der Fehlersuche beauf- 
^ tragt wird^ Mit viel Erfahrung kann man. die Ursachen der Sin- 
10 gularitat auffinden und gegebenenf alls auf Fehler bei der Mb- 
dellierung des technischen Systems zuruckschlieUen. Dies ist 
regelmali'ig *sehr zeitaufwendig und kostspielig. 

■ Es ist daher Aufgabe der Erfindung,- ein Verfahren zuiu Bereit- 
15 stellen von Fehlerinf ormationen tlber Inkonsistenzen in einem 
System von Dif f erentialgleicKungen anzugeben. 

Di'ese Aufgabe wird durch den Gegenstand der unabhSngigen An- 
sprtichG .gel5st . Vorteilhafte Weiterbildungen ergeben sich aus 
20 den jeweiligen Unterarispruchen . ■ 



• Die Erfindung beruht auf dem Gedanken, bereits vor einer- Si- 
■ mulation, . d.h, vor deia AusfUiiren numerischer Berechnungs- 
schritte oder der Zusaramenstellung eines Analogrechners, das 
2 5 System von Dif f erentialgleichungen" mit kombinatorischen Me- 
' thoden auf Inkonsistenzen hin zu untersuchen, Beim Vorliegen 

von Informationen tiber solche Inkonsistenzen kann auf Fehler 
bei der Modellierung des technischen Systems oder auf Fehler 
im technischen System selbst zurUckgeschlossen warden, die in 
30 einem spateren Sxmulationsschritt unerwUnschte AbbrUche 'oder 
gar die Nicht simul ierbarke it des • Systems yerursachen konnen. 
Ein zeitraubendes AusfUhren der Simulation bis zu einem Ab- 
bruch bzw. ein Simulationsversuch ist nicht notwendig. Die 
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. ■ . 3 . 

Korrektur des Modells- des technischen Systems oder des tech- 
nischen Systems selbst wird durch das Vorliegen von Informa- 
tionen uber Inkonsistenzen wesentlich erleichtert . 

5 Das erf indungsgemalie Verfahren sieht in einem ersten Schritt 
vor, eine Abhangigkeitsmatrix A aufzustellen. Die Abhangig- 
keitsmatrix A hat genau so vieie Spalten vyie die Dimension 
des Losungsvektors x des Systems von Dif f erentialgleichungen. 
Die Abhangigkeitsmatrix ^ .hat genau so viele Zeilen wie" die 
10 . Anzahl der. Dif f erentialgleichungen in dem gegebenen System 
von Dif f erentialgleichungen. 



Die Erfindung ist dabei nicht . nur auf Dif f erentialgleichungs- 
systeme der For^ / (^ x ^(0 . > ^''\(^V^)= Q • beschr^nkt . Es kann 
. 15 auch fur Sbnderfalle angewendet werden, bei denen das Glei- 
chungs system zur Simulation des technischen Systems bei- 
' spielsweise die Form f(t,x(t))=0 aufweist. In einem solchen 
-Fall spricht man nicht von einem System von Dif f erentialglei- 
chungen, sondern von einem im .Allgemeinen nicht-linearen 
20 'Gleichungssystem* " ' 

AuJierdem muft nicht unbedingt das zur • Simulation Vorgesehene 
Systeia von Dif f erentialgleichungen selbst auf Inkonsistenzen 
untersucht werden,- Vielmehr kann ein anderes System von Dif- 
^ 25 ferentialgleichungen, welches dasselbe technische System oder 
denseiben technischen Prozefi beschreibt , auf Inkonsistenzen 
untersucht werden. Dies kann vorteilhaft sein, wenn das 
let z^tgenannte .System' von. Dif f erentialgleichungen einen einfa- 
cheiren Aufbau hat als dasjenige, welches spater simuliert 
30 werden soil. Zum Beispxel sind bei der Simulation elektri- 

scher Netzwerke ublicherweise die Gleichungen der modifizxer- 
ten Knotenspannungsanlayse Grundlage der numerischen Simula- 
tion, obwohl die sogenannten Zweigspannungs-Zweigstrom- 
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Gleichungen einen ^esenflich einf acheren Aufbau haben. Durch 
die Unt^rsuchung von Systexae^. von Dif ferentialgleichungen m.t 
einfachem.Aufbau konnen sogar oftmals mehr Inkonsistenzen er- 
mittelt werdeB Oder Inkonsistenzen starker eingegrenzt wer- 



den . 



10 

\ 



Die Abhangigkeitsmatrix A weist dabei Elemente A^j) auf, 
die entweder auf den Wert "Null" oder auf einen Wert "un- 
gleich Null" gesetzt werden, wobei es auf die absoluten Werte 
der Elemente von A nicht ankoinxnt. Dabei wird das Element 
Ait/j) auf einen Wert *0 gesetzt, wenn die i-te- Zeile von f 
vom' j-ten Element von x oder von einer der Ableitungen des j- 
ten Elements von x abhangt. in' alien Ubrigen Fallen wird das 
Element A{i,j). auf den Wert "0" gesetzt. ; , 

Ftlr das erfindungsgemaiie Verfahren ist es unerheblich, ob die 
Abhangigkeitsraatrix A tatsSchlich explicit aufgestellt wird. 
Vieliaehr kommt es darauf an, dafi die entsprechende Informati- 
on, Ob die i-te zeile von f vom i-ten Element von x oder von 
20 einer seiner Ableitungen abhangt, oder nicht, in den folgen- 
den Verfahrensschritten verfugbar ist. 

' • im zweiten Schritt des erf indungsgemaJien Verfahrens wird eine 
Menge von' Zeilenmaschen ' bestiramt, die jeweils die Nummern 
> 25 " derjenigen Zeilen der Abhangigkeitsmatrix A bzw. derjeiiigen , 

Zeilen des Systems von Dif ferentialgleichungen aufweisen, die 
' voheinander abhangig .sind. AuBerdem wird eine Menge von Spal- 
. ^ ' tenmaschen der Abhangigkeitsmatrix a ermittelt, die jeweils 

. die Nummern derjenigen Spalten der Abhangigkeitsmatrix A bzw. 
30 derjenigen Komponenten von x Oder einer seiner .Ableitungen 
aufweisen, die voneinander abhangig sind. 
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Unter einer "Zeileniuasche" wird hier folgendes verstanden. 
Eine Menge Cz nattirlicher Zahlen i, l<><n, heifit. Zeilenma- 
s Che einer Mat.rix mit n Zeilen und m Spalten, wenn sie die 
folgenden Bedingungen erfilllt: 



10 



1-5 



(i) Es gibt keine Transversale T der Matrix A derart, dafi Cz 
in der Menge der Zeilenindizes , von' T enthalten ist. 

(ii) zu jedem Element c von Cz gibt es eine Transversale T . 
von A, so .daB. Cz\{c} ganz in der Menge der Zeilenindizes 
von T enthalten ist. 

Dabei stent der Ausdruck Ca\{c} diejenige Menge dar, die 
erit stent, wenn man aus der Menge Cz das Element c entfernt.- 

Unter einer Transversalen T- der Matrix A wird folgendes ver- 
standen. Eine Transversale. einer Matrix h mit n Zeilen und m 
.Spalten ist eine von gegebenenf alls mehreren moglichen Mengen 
yon Positionen 0, j) von von Null versehiedenen MatrixeintrS- 
gen 4(/V/)..von denen keine zwei oder mehr in der seltaeri Zei- 
le Oder Spalte' stehen. Eine Menge T von Paaren {i>j), 
l</£»,l<y^wi/. stellt eine Transversale der Matrix a dar, 
falls. T die folgenden Bedingungen erfUllt: 

25 (i) fiir alle Elemente von T gilf Aii,j) ist verschieden 

von 0, 

(ii) es seien j) und,(i',/) zwei verschiedene Elemente von ' 
das heiat oder J ^ f r dann gilt und zugleich 



20 



30 



Fur die Erfindung ist es dabei weniger von Bedeutung, bb man 
eine der meglichen Transversalen . der Abhangigkei.tsmatrix A 
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ausrechnet. Die Definition des Begriffs "Transversale" wird 
hier vielmehr benotigt, um die Begriffe "Zeilenmasche" und 
"Spaltenmasche" zu veranschauliQhen . ■ • . 

Unter der Bezeichnung " Zeilenindizes der Transversals T" wird 
folgendes verstanden. Die Menge Z der "Zeilenindizes der 
Transversale T" hat als Elemente die Zeilenindizes der. Ele- 
mente-von T, Mit anderen Worten bedeutet dies, dafi die Menge 
2 genau diejenigen Zahlen i enthalt/- fUr-die es ein" j gibt, 
so daJi das Element (/. v) ein Elemeni: von T ist. 

unter einer Spaltenmasche. wird eine Menge Cs natUrlicher Zah- 
len i, l£/<»», einer Matrix ^ mit n Zeilen und m Spalten 
■ verstanden, wenn sie folgende Bedingungen erftlllt: 

V 

(i) . ES gibt keine Transversale T der Matrix A derart, dafi Cs 

in der Menge der Spaltenindizes von T enfchalten ist, 

(ii) Zu jedem Element c von Cs gibt es sine . Transversale T 
der- Matrix A, so dafl Cs\{c} ganz in der .Menge der Spal- 
tenindizes von T enthalten ist. 



/ Unter der "Menge der Spaltenindizes der Transversale T". wird 

eine Menge Z verstanden, die als Elemente die Spaltenindizes ■ 
der Elemente der Transversalen T der Matrix a enthait- Mit 
,25 anderen Worten bedeutet dies, dafi die, Menge Z genau diejeni- 
# gen Zahlen j enthalt, fiir die es ein i gibt, so. daJJ das Ele- 

■ . ment {i, j) ein Eleiaent der Transversale T der Matrix A ist. 

Nach dem Durchftihren des erf indungsgemafien Schritts 2 liegt, 
30 eine Menge. von Zeilenmaschen vor, die jeweils die Nummern 

derjenigen Zeilen der AbhSngigkeitsmatrix A bzw. derjenigen 
Zeilen des Systems von Dif f erentialgleichungen aufweisen, die 
voneinander abhangig sind- AuBerdem liegt. eine Menge von 



Datum 28.05.01 16:14 FAXG3 Nr: 384577 von NVS:FAXG3.I0.01 01/00 (Seite 8 von 25) 



2S-MAI-2001 IS: 18 PATENTANWALT +49 89 34Q1899B S.09:--35 

■ FIN 101 P/200014448 



Spaltenmaschen d.r I^hangigkeitsnvatrix ^ vor, , die jewexls dxe 
Nuinmern derjenigen Spalten der Abh^ngigkeitsmatrix ^ bzw. 
derj^nigen Komponentan von x oder einer- seiner Able itungen . 
aufweisen, die voneinander abhangig sind. 



. Mit dieser Infornvation kann besonde.s ■ einfach auf Struktur- . 

fehler des als Ausgangspunkt vorhandeneri Dif f erentialglex- 
• chungssyste.s g.schlossen werden, Solche Strukturf ehX.. bzw.. 
xnkonsistenzen in de. urspr^glich.n Dif f .rent.algl.xchungs- 
\ .10 svste. sind h.ufig Ursachen fUr -Fehler bei der.Berechnung de. 
■ Losung des Differentialgleichungssystems. 

Ein grundlegender Schritt beixa Ausfiil^ren des erf indungsgem^- ' 
' fien Verfahrens liegt in dem Auffinden von Zeilemnaschen und 

15 spaltenmaschen. Verfahren zur. Auffinden von Zeilenraasch.n und 
Spaltenxaaschen sind bekannt- In; dar Literatur finden -sich un- 
ter dem Stichwort "minimally structurally singular subsets of 
equations" (vgl. CC Pantelides, The consistent in.tn.alxza- 
tion of ■ differential-algebraic systems in SIAM J. Sci. Sta 
20 tist. Comput./ 9(2):213-231,,Marz 1988), Verfahren- zur Be- 
re'-hnung von " Zeilenmaschen" . Zum Bestiimen von "Spaltenma- 
schen" kann auf ein Verfahren zum Bestixnmen von "Zsilenma- . 
schen" zurllck gegriffen werden, werin von derjenigen Matrix, 
von der die "Spaltemuaschen" bestiitant werden sollen, eine 
25 transponierte Matrix bestiirout wird. Zum Bestimmen der Spal- 
'^^ tenmaschen der Matrix werden dann die Zeilenma schen , der- ent- 

~ ' sprechenden trahsponierten Matrix bestimmt. , 

• Diese Art der Bestiininung der Zeilenmaschen und ^paltenmaschen 
30 fur das erf indungsgemSiie Verfahren soil nicht beschr^nkend 
• verstanden werden. Vielmehr ktinnen auch andere Verfahren zui. 
Bestimmen von Zeilenmaschen und Spaltenmaschen verwendet wer- 



Datum 28.05.01 16:14 FAXG3 Nr: 384577 von NVS:FAXG3.l0.0101/00 (Seite 9 von 25) 



■28^MftI-200i 16: 18 " PflTENTflNUJALT • +4y 

FIN 101 P/200014448 , 



. ' ■ den. 3pfern die Ubrigen B.dingungen far da. Vorhandens.ein von 
Zeilenmaschen .und Spaltenmaschen erfiiHt werden. 

im abschli^Benden Schritt des erf indungsgemSfien Verf ahrens 
5 warden fUr jede Z^ilenmasche und.f^ir jede Spaltenxaasche, di« 
gema^ Schritt 2 bestimzat worden sind, . di0 darin enthaltenen 
'■ Numrcern ausgegeb^n. Diese Nummern deuten im Fall der Zeilen- 
maschen auf die. laufenden Nummern der qieichung'en des Dxffe- 
. . rentialgleichungssystems bin,, die moglicherweise mit einem 
10 strukturellen -problem behaftet sind. Im Fall der- Spaltenma- 
. schen deuten die darin enthaltenen ^nd als Fehlerinf ormatxo- 
■ hen ausgegebenen Nummern auf die Nummern der Komponenten des 
Losungsvektors x bin. die moglicherwei^e mit einem struktu- 
reilen Probiem bfehaftet sind. 



15 



20 



25 



4 



30 



Mit dem erfindungsgemafien Verfahren ist es mbglich, schon vd.r 
dem Ausfuhreri von numerischen Berechnungsschritten Oder der 
ZusanB^venstellung einas AnaXogrechners die Struktur, de^s zu be- 
rechnenden Systems von Dif f erentialgleichungen auf ■ seine Kon- 
sistenz bin zii Uberprlifen. Wenn mit dem .erf indungsgemSfien 
verfahren ein System von Dif f erentialgleichungen, als inkpnsi- 
stent erkannt worden ist, konnen sine Grofizahl der mSglichen 
struktufell.en Tehler des Systems von Dif f erentialgleichungen 
auf gefunden werden. Dadurch wird die Fehlersuche beschleu- ' 
nigt. Zudem werden zeitraubende Simulationsversuche verrtiiey 
den. • ■ . • . ' ' . 

Das erf indungsgemafie verfahren ist bei jegiicher Art eines 
.Simulationssysteitis zur numerischen LGsung von Differential- 
gleichungssystemen sowie zur Losung von Differentialglei-. 
chungssystemen mittels Analogrechner anwendbar. 
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Bei einer besonders vorteilhaf ten Weiterbildung der Erfindung 
wird vor dem Ausfuhren von Schritt 1 des erf indungsgemaSen 
verfahrens eine Gleichungs-Bedeutungsliste der Lange n ange- 
legt, in der jeder Gleichung des Gleichungssystems eine Glei- 
5 chunasnuimner und/oder eine Gleichungs-Textinf ormation zuge- 
ordn.t wird. Genauso wird vor dem Ausftihren- von Schritt 1 des 
erfindungsgemalien Verfahrens eine Komponenten-Bedeutungsl.s^e • 
' der L^nge m angelegt, in der jeder Koiuponente des LGsungsvek- 
tors X eine Koiaponentennu^mer und/qder eine Komponenten- 
10 Textinformation zugeordnet wird. Dabei wird die.jeweils in 

der Bedeutungsliste abgelegte Gleichungs-Textinf ormation bzw. 
\ . Komponenten-Textinformation vorteilhaf terweise so ausgewahlt. 

dafi diese eine Qedeutung im Zusainmenhang mit dem zu simul-ie- 
. renden technischen System erhalten. Dadurch konnen Bestand- 
15 teilen des' simulierten technischen Systems und Teilstrukturen 
des zu simulierendeh technischen Systems Gleichungen des Dif- 
f erentialgleichungssystems und Komponenten des LOsungsvektors 
■ X zugeordnet warden, die eine Interpretation der mit dem er- 
,, . • . f'indungsgemaaen Verfahren ausgegebenen Fehlerinformationen 
20 erleichtern. • ' 

in schritt 3 des erf indungsgemalieh Verfahrens wird dazu an , 
der Stelle des Ausgebens der in jeder Zeilenmasche enthalte- 
nen Nmnmern das Ausgeben der Gleichungsnummer und/oder der 
25. Gleichungs-Textinformation entsprechend der Gleichungs- 
M • Bedeutungsliste vorgesehen. Fiir den Fall, dafi in der Zeilen-- 
masche eine gan.ze Zahl i enthalten ist, wird nicht die Nummer 
i ansgegeben, sondern vielmehr der Inhalt der i-ten Komponen- 
te der Gleichungs-Bedeutungsliste. Genauso wird in Schritt 3 
30 ■ des'erfindungsgemaiien Verfahrens. an der Stelle des Ausgebens 
der in jeder Spaltenmas.che- enthaltenen Nmnmer die Komponen- • 
tennummer und/oder die Komponenten-Textinf ormation entspre- 
' chend der Komponenten-Bedeutun-gsliste 'ausgegeben. . Fur den- 
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Fall, daU in der ^Spaltenmasche die ganze Zahl j als Nuinmer 
vorhanden ist, wird. nicht die .Nummer j ausgegeben, sondern . 
vielinehr der Inhalt der j-ten Komponente der Komponenten- 
Bedeutungsliste . 

5 ' ■ ■ . 

Durch diese Ausgestaltung des erf indungsgemafien Verfahrens . 
. werden als Fehlerinf ormationen Bedeutungsgehalte ausgegeben, 
die einen unmittelbaren Bezug zum zu simulierenden- tedhni- 
schen system haben. Mit diesen Fehlinf ormationen konnen. sy- 
^ 10 stematische Fehler In der Struktur des Systems von Diff eren- 
tialgleichungen zur Beschreibung des technischen Systems be- 
sonders leicht veranschaulicht werden, so dafl die Fehiersuche 
weiter beschleunigt wird, 

15 ' Bei den erf indungsgemaJien Verfaihren wird ein vorzugsweise di- 
, gitaler .Rechner eingesetzt, der zvamindest einen Speicher, ain 
Rechenwerk, eine Eingabevorrichtung und eihe. Ausgabevorrich- ■ 
tung. aufweist . . . 

.. 20 Die Erfindung ist auch in einem Computerprogramm zum .Bereit- 
■• stellen.von Fehlerinformationen Uber Inkonsistenzeh in einem 
Gleichungssystem verwirklicht. Das Computerprogramm ist dabei 
so ausgebildet, daJJ nach Eingabe der Systemeigenschaften, der. 
Start- Oder Randbedingungen und der Einf lufinahmen auf das Sy- 
25 . stem ein erf indungsgemafies Verfahren gemafi einem der vorher- 
gehenden Anspruche ausfxihrbar ist. Dabei kSnnen Ergebnisse. 
einer Simulation, der L6sungsvektor oder die Losungsvektoren 
zu verschiedenen Zeitpuiikten ausgebbar sein. Es ist aber auch 
mbglich, daB das Computerprogramm lediglich Informationen 
30 iiber Inkohsistenzen gemaii einem der vorhergehenden Ansp ruche 
. bereitstellt .' 
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Durch" das erf indungsgemalie Computerprograiron ergeben sich er- 
hebliche Laufzeitverbesserungen gegenObe. dsn bekanrxten Simu- 
lationsprogrammen und Simulationsverf ahren, well zahlrexche 
f.hlerhafte Programmlauf e bzw. Simulationsversuche vermxed.n 
werden konnen, 

Oie- Erf .ndung betrif ft au.erdexa einen Datent^r.ger .it einem 
sclchen computerprograxnm'.owie ein Verfahren, bei dem- exn 
solch.s Conputerprogra^ aus einexn elektronischen Datennetz, 
wie beispielsweise aus dem Internet auf einen an das Daten- . 
netz angeschlossenen Computer heruntergeladen wxrd. . 



10 



15 



20 



Das erflndungsgemaae/Computersystem ist so ausgebildet, daB . 
auf die.em ein erfindungsgem.Bes Verfahren zum Bereitstellen 
von Fehlerinformationen Ober Inkonsistenzen • in einem- Glex- 
chungssystem. ausfUti'rbar ist, . . 

Die Erfindung betrif ft schlieBlich auch die . Verw.endung aines 
erfindungsgera^I^en Verfahrens und/oder eines Computersystems 
zum Bereitstellen von" Fenlerinformationen uber Inkonsistenzen 
in einem Gleichungssystem. 



25 



Bei dem erf indungsgemafien Verf ahren wird. dabei zun^chst exn 
das system- Oder den'Prozefi beschreibendes Dif f erentialgle^-- 
_ chungssystem oder andere Formen wie beispielsweise nicht- 
^ M. linaare Gleichungssysteme, di4 keine Dif f erentialgleichungen . 

^' aufweisen, aufgestellt. Dabei wird hiar davon ausgegangen, 

dafi dem Fachmann die liierzu notwendige Vorgehensweise bekannt 
ist. in der in der Beschreibungseinleitung genannten Fachlo. ^ 
teratur wird hierauf ausfuhrlich eingegangen. Im iibrigen exi- 
stieren hierzu-einschlagige Computerprogramme. 



30 
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12 

Das erste AusfUhrungsbeispiel gemali der Erfindung betrifft 
die Simulation eines hier nicht geseigten technischen Sy- 
stems, das durch numerisches Losen des Gleichungssystems 



10 



der Form 



' /3 (:<:, (4^. (4 -^, (4 *. (4*.(')t*3(4A.i'../'s) = 0' 
tuit den. Funktionen fl, f2 und f3 und mit dem Parametervektor 
p gemaB Figur 1 beschrieben warden kann. 

Zur vorhersage des verhaltens des Systems wird das Glei- 
15 chungssystem numerisch gel6st, dM. es werden Werte fUr den 
Unbekanntenvektor x(t) zu einem odef mehreren Zeitpunkten t 
berechnet. Hierzu werden hier nicht .veranschaulichte bekannte 

■ numerische Losungsyerf ahren eingesetzt, die auf einem hier 

■ nicht gezeigten Computersystem als Computerprogramm ablaufen. 

20 ' . _ 

• ErfindungsgemaB werden vor dem eigentlichen Losen des Glei- 
chunssystems Schritte zum Bereitstellen- yon Fehlerinf ormatio- 
nen uber Xnkonsistenzen im Gleichungssystem ausgefUhrt. 

25 Im Schfitt 1 des erfindungsgemafien Verfahrehs wird die in Fi- 
gur 2 angegebene Abhangigkeitsmatrix- A bestimmt, wobei alle 
von Null yerschiedenen Elemente der Abhangigkeitsmatrix h in 
Figur 2 durch einen Stern ("*") gekennzeichnet sind. 
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Das beispielhaft herausgegrif f ene Element. ^ (l, 1 ) wird auf 
einen im wesentlichen beliebigen- von 0 verschiedenen Wert 
gesetzt,. denn das erste Element von t, also 

5 A(t,x{t),x(t),...,x"''it)) , hangt vom ersten Element von x, also von 
' •. xi (t) ab. ■ ' 

Das beispielhaft herausgegriff ene Element 4(3,1 ). wird auf. 
den Wert "0" gesetzt, denn das dritte . Element von f, also 
10 (0) r i3t vom ersten Element von X, also von 

■xi(t) und von den Ableitungen des. ersten Elements von. X, also 
''kk von xi'^'(t), unabhangig. 



15 



20 



im schritt 2 des- erf indungsgemafien Verfahrens wird mit dem 
Ziel, zeilenmaschen und Spaltenmaschen der in Figur 2 angege- 
beWen Abhangigkeitsmatrix E zu ermitteln, zunachst die in Fi- 
gur 2 angegebene Transversals. T der Abhangigkeitkmatrix h be- 
stimmt. Verfahren zum Ermitteln von Transversalen (von "maxi- 
mum ' cardinality bipartite matchings") sind in 1... S. Duff, "On 
algorithms for obtaining a maximum, transversal" , ACM Trans. 
Math, software^ 7 (3) : 315-330, 1981, in E . L. Lawler, "Combi- . 
natorial Optimization: Networks and Matroids",. Holt, Rinehart- 
and Winston., 1976, in L. Lovasz and M. D. Pl{immer, "Matching 
Theory", North-Holland Mathematics Studies 121. Annals of 
25 Discrete Mathematics, 29. North-Holland," 198 6, Oder in C. C. 
Pantelides, '"The consistent initialization of differential- 
algebraic Systems", SIAM J. Sci. Statist. Comput . , 9 (2) : 213- 
231, Mar. 1988, angegeben. ' 



30 



Die Transversals T aus Figur '2 hat 3 Elements. Da n - m - ,3 
ist," gibt es weder Zeilenmaschen, noch Spaltenmaschen der Ab- 
hangigkeitsmatrix A a.us Figur 2. Das in C. C. Pantelides, 
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10 



15 



20' 



25 



30 



"the consistent initialization of differential-algebraic Sy- 
stems", angegebene Verfahren zuiu Berechnen von Zeilemuaschen 
("minimally structurally singular subsets of equations") fxn- 
det keine solche Zeilenmasche der Abhangigkeitsmatrix A.. Das 
dort angegebene Verfahren zum Berechnen vo.n zeilenmaschen 
findet auch keine solche Zeilenmasche der Transponisrten der 
Abh^ngigkeitsxuatriK A. Folglich gibt es keine Spaltenmasche- 
de^ Abhangigkeitsmatrix A. Die in Schrit.t 2 besti^ute Menge 
von zeilenmaschen .der Abhangigkeitsmatrix ^ ist leer. Die m 
--hritt 2 bestimmte Menge von Spaltenmaschen der AbhSngig- , 
keitsmatrix h ist • ebenf alls leer. I.m Schritt ^3 des erfin- 
dungsgemaiien .Verf ahrens wird demnach keine Fehlerinformatxon 
ausgegeben. 

Das eingangs angegebene Gleichungssystem kann mit hoher Wahr^ 
scheinlichkeit. ohne das Auftreten von Fehlern gelost . werden , 
Wird die Simulation dennoch abgebrochen, oder liefert sie 
nicht plausible Losungen, oder ist sie ganz unmoglich, dann 
konnen strukturelle Fehler in der Beschreibung des. techni- . ■ 
schen systems oder'im technischen system selbst mit hoher 
wahrscheinlichkeit ausgeschlossen werden. Dies vereinfacht 

die Fehlersuche erheblicti, 

t ' * ' ■ 

Das zweite AusfUhrungsbeispiel gemafi der Erfindung betrif.ft 
die Simulation eines weiteren hier nidht' gezeigten techni- 
' schen systems, dessen Verhalten. durch numerisches Losen des 
Gl e i Chungs s ys terns 

der Form ■ , 
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mit den Funktionen xl, f2 und f3 und irit. dem Parametervektor 
, £ gemaib Figur 3 beschrieben werden kann-- 

Zur Vorhersage des Verhaltens des Systems soli das Glei- 
chnngssysteiu numerisch gel6st warden, d.h. es soll.n Werte 
fUr den Unbekanntenvektor x(t) zu ein^m oder mehreren Zext- 
punkteht b^rechnet werden. Hierzu werden hier nicht veran^ 
10 schauliehte b^kannte n^erische L5sungsverf ahren eingesetzt, 
die auf eineiu hier nicht gezeigten Computersystem als Gompu- 
t.erprdgramm abiaufen. ' 

■ Erf indungsgemaB werden vor dem eigentllchen. L5sen des Glei- 
.5 chunssystems ■ Schritte zum Bereit.tellen von Fehlerinformatxo- 
• nen ilber Inkonsistenzen im Gleichungssystem ausgefUhrt. 

im Schritt 1 des erf indungsgemafien Verfahrens. wird die in Fi- 
gur 4 angegebene Abhangigkeitsmatrix & bestimmt, wobei alle 
20 von Null verschiedenen Elemenfe der JUDhangigkeitsmatrix h ^n. 
Figur 4 durch einen Stem ("*") gekennzeichnet sind. 

Das beispielhaft herausgegriffene Eleiaent 4 (.1,1.) ^ird auf 
. einen im wesentlichen beliebigen von 0 verschiedenen Wert. 
" 25 gesetzt, denn die erste Zeile von f , also 

^ ^<^')' - ' ' ■ ersten Element von x, also von 

xi(t) "ab. Dasselbe gilt fUr die Elements ^4(1, 2)-, di^^) > 
Ai2,3 ) und v4(3,3 ) -. 



Datum 28.05.01 16:14 FAXG3 Nr: 384577 von NVS:FAXG3.I0.01 01/00 (Seite 17 von 25) 



28-MA 1 -2001 • 16=21 PflTENTflNWfiLT . +49 aH Jrfmayyvd 

FIN 101 P/20001444a 

16 • • 

Das beispielhaft herausgegrif f ene Element . 4 ( 3, 1 ) wird auf 
den Wert "0".gesetzt, denn das dritte Element von f, also ' .. 
y3tx(^U(/),-.i'*'(0) . vom ersten Element von x, also von 

xi(t) und von den Ableitungen des .ersten Elements von' x, also 
5 von xx'^'(t), unabhangig. Dasselbe gilt ftir die Elements.. 
4(2,1), 4(3,2). und 4(2. 2.) . 

im Scliritt 2 des erf indungsgemaBen Verfahrens .wird mit dem 
• Ziel, Zeilenmaschen und Spaltenmaschen der in Figur 4 angege- 
10 benen Abhangigkeitsmatrix A ermitteln, zunachst die in Fi- 
gur 4. angegebene Transversals T der -Abhangigkeitsmatrix A be- 
stimmt. 

t. . • " , . ■ - , 

Das in C. C. Pantelides, "The consistent initialization of 
15 differential-algebraic Systems", angegebene Verfahren zum Be- 
rechnen von Zeilenmaschen ("minimally structurally singular 
subsets of : equations") findet die Zeilenmasche {2. 3} der Ab- 
hangigkeitsmatrix Die Menge Z der in Schritt 2 gefundenen 
Zeilenmaschen der Abhangigkeitsmatrix h ist in Figur 4 ange- 
20 geben. 

* ' . ■ . . ' 

Das in e. C. Pantelides, . "The consistent initialization of 
differential-algebraic Systems", angegebene Verfahren zum Be- 
rechnen vdn Zeilenmaschen a:ngewendet auf die Transponie.rte 

25 der Abhangigkeitsmatrix ^ findet die Zeilenmasche {1, 2 } der 
■ Transponierten der Abhangigkeitsmatrix a: Diese Zeilenma- 
. sche der Transponierten der Abhangigkeitsmatrix ^ wird erfin- 
dungsgemaa' als Spaltenmasche der Abhangigkeitsmatrix A aufge- , 
•faBt. Die Menge S der in Schritt 2 gefundenen Spaltenmaschen 

30 der Abhangigkeitsmatrix ^ ist in Figur 4 angegeben. Im 

Schritt 3 wird die in Figur 5 angegebene Fehlarinf ormation 
ausgegeben.- . ' • 
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GemaB der Erf indung .wxrd eine Vorhersage des Verhaltens • des 
zugrundeliegenden System, mit nuxaerischen Mitteln gar nicht 
erst versucht, well dabei Fahler auftreten werden, wenn exne 
Simulation uberhaupt moglich ist. Vielmehr. mUssen die Model- 
lierung das Systems und das System selbst noch einmal uber- 
pruft werden. Dadurch wird' wertvollB Rechenzeit auf dem hxer 
nicht gezeigten Cpmputersystem eingespart- Die Fehlersuche ' 
wird durch die Kenntnis der in Schritt 3 ansgegebenen Fehle- 
rinformationen erheblich vereinf acht . 



I'O 



15 



Ein drittes Ausf uhrungsbeispiel gemafi der Erfindung betrifft 
das in Fig-. 6 gezeigte technische System, dessen Verhalten 
durch numerisches Losen des Gleichungssystems 



der Form 



'20 



/, (4 (4 (4 a./'2.a)= 0 . 

mit den Funktionen fl, f2 und f3 und mit dem Parametervektor 
£ gemaJi Figur 3 beschrieben werden kann. 

Zur Vqrhersage des Verhaltens des Systems soli das Glei- 
25 chungssystem nmtverisch ge lost werden, d.h. es sollen Werte 
fair den Unbekanntenvektor x(t) zu einem oder mehreren Zeit- 
punkten t berechnet werden. Hierzu werden hier' nicht veran- 
schaulichte bekannte n;imerische Lbsungsverf ahren eingesetzt, 
. ''. die auf einem hier nictit gezeigten computersystefti als Compu- 

'30 terprogramm ablaufen. . 
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Erf indungsgex^a. warden dabe. vor dem eigentlichen L5sen des ; 
■ ^^^ys.^s Pe.ler...or.ationen ..e. Xn.onsiste.z.n 
~ betreffend.n Gleichungssystem wie folgt be.ei.tgestellt.. 

5 Die Lasungen des betreff.nderx Qleichungssystems sind Ruhela- 

g.n bzw. "Arbeitspunkte'- oder. "DC-L5sungen" des in Fxgur 6 - 
. .angegaWen ele.trischen KetzwerKs,. w.lche^ aus den folgenden 
Netzwerkelementen ( "Bauelementen" ) besteht: . 

- einem iinearen Wid^rstand mit dem Widerstandswert R zwx- 
10 schen den Knoten 1 und 2 des Netzwerks, 

- einer linear^n Kapazitat mit dem Kapazitatswert CL zw.- 
schen den Knoten 1 und 0 des Netzwerks, und 

. - einer linearen Kapazitat mit dem Kapazitatswert C2 zwx- 
spben den Knoten 2 und 0 des Netzwerks. 

Daraus ergeben sich die Koraponenten de. in Figur 3 angegabe- 
nen Parametervektors p entsprechand den Werten CI, 1/R,- C 2 . 

Dxe Komponanten xl(t), x2(t) und x3 ft) voH x(t) entsprechen- 
20 den folgenden Graven des in Figur 6 angegebenen .Netzwerks: 

- xl-(t) entspricht der Spannung zwischen den' Knoten 1 und 

, 0, ' '. ■ • . • . 

- x2(t) entspricht der Spannung zwischen den Knoten 2 und 

0, und 

25 - x3(t) entspricht der Spannung zwischen den Knoten 1 und 
'2'. .. • • . ... ■ 

' Die erste Gleichung des Gleichungssysteius, also die Gleichung 



30 



fl(xl(t), x2(t), x3(t), pl, P2, p3) - 0 
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, Bi, zweite .leichung des Gleiohungs.yste»s, also die 01=1- 

Chung - . " . " 



.f2(xl(1i), ^2(t), X3(t), pl, P2. P3) - 0 

10 '' 



xs 



t die Kirchhoffsche Stron>gleichung • f tar den Knoten 1 des 



, .iSTetzwerkes. aus Figur 6. 



•15 Die dritte Gleichung des Gleichungssystems, also die Glei-, 



Chung 



20 



•f3(xl{t) , x2(t) , .x3(t), pl, p3) - 0. 

ist die Kirchhoffsche Stromgle.ichung fur den Knoten 2 des^ 
Netzwerkes aus Figur 6. . 

. ~ E.«.duna.,em.s= wl.d die In Fi^r 7 .„,.=,eb.ne Glei=hu„,.b.- 

dputun,3li3« G u«d die in figur 7 angeg^ben. Ko«pon.ntenba- 
deutungsliste K angelegt. 

Wie in den vorhergehehden Ausfiihrungsbeispielen wird die in 
30 Figur-4 angegeben. Menge Z der eriaittelten Zeilenma^chen and 
die in Figur 4 angeg.bene Menge S der ermit-telten Spaltenma. 
■ schen bestimmt. 
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im Schritt 3 des erf indungsgemaJJ weitergebildeten Verfahrens 
wird unter Vsrwendung der Gleichungsbedeutungsliste G und der 
Komponentenbedeutungsliste K aus Figur 7 die in Figur 8 ange- 
gebene Fehlerinformation ausgegeben. 

AUS der in Figur 8 angegebenen Fehlerinformation laiit si.ch 
erkennen. dafidie Kirchhof f schen Stromgleichungen zu den Knp- 
ten -r und 2 des Netzwerks aus Figur 6 voneinander linear ab- . 
hangig sind und daB di^ beiden Spannungen zwischen Knoten 1 
und 0 und zwischen Knoten 2 und 0 des Netzwerks aus Figur 6 
fiir Ruhelagen dieses Netzwerks nicht eindeutig bestiiorat wer-. 
. den k:5nnen- 



10 



15 



20 



GemaB der Erfindung wird eine Vorhersage des Verhaltens des- 
zugrundeliegenden Systems mit numerischen Mitteln nicht vor- • • 
genoitaaeii, weil dabei- Fehler auftretcn werden, wenn eine. Sxmu- 
•lation tiberhaupt mdglich ist. vielmehr ■ mussen die Modellie- 
rung des Systems und das System selbst noch einmal tlberprUft. 
werden. Dadurch wird wertvolle Rechenzeit auf dem hier nicht. 
gezeigten Computersystem eingespart . Die Fehlersuche wird 
durch die Kenntnis der im Schritt 3 ausgegebenen Fehlerinfor- 
mation erheblich vereinfacht. 
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Patentanspruche • 

1 V.rf ahren zum Beraitstellen von Fehlerinf ormationan ub«r 
■ inl^onsistenzen in einem mit .Lne^ auf einezn Computer ab- 
■ laufenden Computerprograrom oder mit einexa Analogrechner 
zu losenden Gleichungssystem der Form - . 



10 



. wobei x(t) und seine Bi=leitungsn x (0,-.:, x^'^ (Oj swells m 
' Elemente aufweisen und wobei p ein Parametervektor, ist, 

• der im Gleichungssystem vorJcommen kann, 

wobei das Verfahren die folgenden Schritte aufweist: 

15 ■ Schritt 1 : , •' 

Aufstellen einer Abhangigkeitsiaatrix lait m Spalt^n und 

mit n Zeilen, und zwar so, dafi 

das Element A{i:j)^0 gesetzf wird, wenn die i-te Zeile 

von f, also 

20 a) vom j -ten Element , von X/ also von (t) , o.der 

b) von einer der Ableitungen des j-ten Element von' x,.. al- 
so von x^'" (t) , abhSngt, ' ; 
und dali sonst das Element" ^(/, 7 )= 0 gesetzt wird. ■ 

schritt 2: 

25 ■ . Ermitteln einer Menge von Zeileniaaschen, die jeweils die 
. Nummern derjenigen Zeilen der Abhangigkeitsmatrix a auf 

• weisen, die voneinander. abhangig sind, und einer Menge 
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von spaitenxaaschen der Atohangigkeitsmatrix a, die jeweils 
die Nummern derjenigen Spalten der Abhangigkeitsinatrxx A 
aufweisen/die voneinander ' abhangig sind/ sofern solche 



Zeilenmaschen bzw. Spaltenmaschen vorhanden sind, . 

5 Schritt 3; . . 

Ausgeben von Fehlerinf ormationen, und zwar fUr jede xm 
■ schritt 2 ermittelte Zeilenmasche und flir jede im Schritt 
2 ermittelte Spaltenmasche die darin enthaltenen Nuironern, 

10 2. Verfahren nach j^spruch 1, 
dadurch gekennzeichnet , dali 

vor dera AusfUhren von- Schritt 1 eine Gleichungs - 
Bedeutungslisteder Lange n.angelegt wird, in der D^der 
". Gleichung des Gleichungssystems eine Gleichungsnummer 
15 und/oder eine Gleichungs-Textinf ormat ion zngeordnet wird, 

• und daa vor deiu Ausfuhren von Schritt 1 eine Komponenten- 
■ Bedeutungsliste der L^nge m angelegt wird, in der jeder. 
. ■ . Komponente des LOsungsvektors x ^ihe Komponentennrnnmer 
und/oder eine Komponenten-Textinformation zugeordnet 

20 - wird/ 

wobei in schritt 3 an Stelle des Ausgebens der in jeder 
Zeilenmasche enthaltenen Nummern das Ausgeben der Glei- 
chungsnmniuer und/oder der Gleichungs-Textinf ormation ent- 

■ sprechend der Gleichungs-Bedeutungsliste vorgesehen ist, 

25 und wobei in Schritt 3 an Stelle des Ausgebens der in je- 

' der Spaltenmasche enthaltenen Nmmuern das Ausgeben der 
• Komponentennuinmer und/oder der Komponenten Textinf ormati- 
■on entsprechend 'der Komponenten-Bedeutungsliste vorgese- 

hen' ist . " 

'3.0 . 

3.. computerprograitanprodukt sowie Computerprogramn zum Be- 
reitstellen von Fehlerinf ormationen tibe'r Inkonsistenzen 

■ • in einera mit einem auf . einem Computer -ablaufenden Coiupu- 
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terprograinia oder mit einera Analogrechnar zu losenden 
Gl.ichungssystem, das so ausgebildet ist, dafi ein Verfah- 
ren gemSli einem der vorhergehenden i^spriiche ausfiihrbar 
ist, * ' . 

' 4. Datentrager .it einem Computerprogranuuprodulct bzw.- Cozapu- 
terprogramm nach Anspruch. 3 . 

5 Verfahren, bei dem ein Computerprograininprodukt b^w. Com- 
10 ■ p.uterprograinzn nach Anspruch 3 aus einem elektronischen •, 

Dat«nnetz wxe beispielsweise aus dem .Internet auf exnen 
an dks Datennetz angeschlossenen- computer beruntergeladen ■ 



15 6, 



20 



7 . 



wird. 



Computersystera, auf dem. ein Verfahren zum Bereitstellen 
von Fehlerinformationen uber Inkonsistenzen in exnem. 
Gleichungssystem nach Anspruch 1 oder 2 ausfUhrbar ist.. 

Verwendung eines Verf ahr ens ' nach ^Anspruch 1 oder 2 
und/ode,r eines. Computersystems nach Anspruch 6 zur Vor- 

tr=»^T, = n-on« Pines Svstems .ausqehend von vorge- 
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Zusawmenf assung 



Verfahren zuiti Bereitstellsn von Fehlerinf ormationea Otoer In ' 
konsistenzen in einem System von pifferentialgleichungen . 

Die Erfindung betrifft ein Vsrfahren. zum Bereitstellen von 
Fehlerinformationen uber Inkonsistenzen in einem- ein techni 
sches syslieiu Oder einen technischen ProzeB beschreibenden Sy- 
stem von Dif ferentialgleichungen. . 

[Fig. 6] , . 
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Fig - 1 ■ 
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* 0 * 
V 0 * * . 



- {(1, 1), (2, 3), (3, 2)} 



Fig * 2 
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A (X, (t), X, (t), X, it), P,,P.,P.>- P2 ■ ^3 



p = 



'p? 

P2 
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0 0 * 



{ a/ 1) ' (2, 3) ) 



Z = { (2, 3}.} 



S - {{1, 2}} 



Fig. .4 
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Eine Zeilenmasche mit den Elementen 
2,3 ■ 

vjurde cfefunden. 

Eine spaltenmasche mit den Eleiaenten 

1/2 . . ^ 
wurde gefxinden. 
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.-l-Kirchhoffsche Spannungsgleichung ftir die aus Cl, 

...ii^i^t^^^^^ den Knoten ... 

-rrKirchhoffsche Stromgleichung ftir den Knoten 2 ) 



K=( "l:Spannung zwischen Knoten 1 und 0^, 
'•2;Spannung zwischen Knoten 2 und 0 , 
"3:Spannung zwischen Knoten 1 und ^ ) 



Fig - 
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"2:Kirchhoffsche S.tromgle^chung f^r a 
"3:Kirchhoffsche stromgleichung fur ae 
wurde gefunden- vip,menten 

"IrSpannung swischen Knoten 1 und 0__, 
"2:spannung.zwischen Knoten- 2 und 0 . 
wurde gefunden. > • 



Fig. 8 
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